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We present a theoretical investigation into the two-center interference in aligned HJ. The influ- 
ence of the laser field on the recombination step is investigated by comparing laser- induced harmonic 
generation with harmonic generation from field-free collisions of Gaussian wave packets with the 
core. We find that for different Gaussian wave packets colliding with the molecule, the interference 
minimum occurs at the same alignment angle. The same result is obtained for the laser-induced 
spectrum when only a single electronic trajectory per harmonic contributes. When multiple elec- 
tronic trajectories contribute, we find an effect on the minimum position because the interference 
between short and long trajectories is alignment-dependent. The two-center interference and the 
influence of the Coulombic potential are clearly seen not only in the harmonic intensity and phase 
but also in the polarization direction and ellipticity. We observe signiflcant ellipticity of the emitted 
radiation around the two-center interference minimum. 

PACS numbers: 33.80.Rv, 42.65.Ky 



I. INTRODUCTION 

When a gas of atoms or molecules is subjected to a 
stronglaser field, high-harmonic generation (HHG) takes 
place [l| . This process converts many of the laser photons 
into a single high-frequency photon in the extreme ultra- 
violet (XUV) or soft x-ray regime. The generation pro- 
cess can be understood in terms of the three-step model 
0. First, the electron tunnels out under the influence 
of the laser field, and then it propagates freely in the 
laser field and can be driven back to the nucleus, where 
it can finally recombine. HHG has been used to gen- 
erate coherent XUV radiation 0, 0] . Concerning HHG 
from molecules, there has been a lot of attention of the 
community recently toward the relationship between the 
emitted radiation and molecular characteristics. One can 
for instance determine the internuclear distance in the di- 
atomic molecules H J , H2 and to a lesser extent CO2 and 
O2 from a destructive interference minimum in the har- 
monic spectrum 0-0| or reconstruct molecular orbitals 
from the emitted radiation d, Q . The correct explana- 
tion of the minimum in CO2, however, appears to require 
a multiorbital treatment [lO[ . Recent experiments on N2 
indicate multiorbital contributions as well Prepar- 
ing molecules to control the harmonic emission is also 
possible: ring current states were used theoretically to 
generate circularly polarized high harmonics [l2j . 

In the three-step model or its quantum-mechanical 
formulation in terms of the strong-field approximation 
(SFA) [31 J due to wave-packet spreading the on-axis con- 
tinuum wave packet carries practically no signature of the 
ground-state wave function at the moment of recombina- 
tion. Therefore, the molecular characteristics imprinted 
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on the emitted spectra must come from the recombina- 
tion step (or from multiorbital interference) . In the exist- 
ing practical molecular imaging techniques the influence 
of the laser field on the recombination step and often also 
the effect of the Coulomb potential on the propagation 
step arc ignored. To study the effects of these approx- 
imations, we perform a numerical comparison between 
the harmonics emitted in a normal laser-induced HHG 
process and harmonics emitted when an artificially pre- 
pared wave packet collides with the molecular ion in the 
absence of any laser pulse Additionally, we compare 
the harmonics generated by both physical and artificial 
laser pulses to disentangle the influence of the different 
harmonic trajectories. We show that elliptically polar- 
ized HHG radiation from linearly polarized generating 
pulses occurs near the two-center interference. Elliptical 
polarization of harmonics from linearly polarized gener- 
ating pulses was also found in recent experiments [14| . 



II. METHOD 

In this article we solve the time-dependent Schrodinger 
equation (TDSE) numerically in two dimensions for a 
molecular ion with a single electron. We consider two di- 
mensions because in 2D many of the 3D characteristics of 
HHG are already present, such as the existence of direc- 
tions perpendicular and parallel to the laser polarization 
axis. On the other hand, the TDSE can be solved very 
quickly, allowing the TDSE to be solved for many align- 
ment angles of the molecule in the laser fleld. The TDSE 
is solved using the split-operator method [H, [l^ , and 
the ground-state wave function is found by imaginary- 
time propagation |17|. 

Atomic units are used throughout this article unless 
indicated otherwise. We will focus on 2D H^ with flxed 
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nuclei. The TDSE reads 

^'-%^^mr,t), (la) 

H='^ + V{r)+r--E{t), (lb) 




(Ic) 

where p is the momentum operator, r = , y is a soft- 
core potential, and E(t) is the time-dependent electric 
field of the laser pulse. The internuclear axis R makes an 
angle 9 with respect to the laser polarization axis x. Us- 
ing the softcore parameter = 0.5 and an internuclear 
distance of i? = 2, the ionization potential is Ip = 30.2 
eV. The total harmonic emission spectrum including x- 
and 2/-polarization is calculated from the numerical solu- 
tion of the TDSE as 

S'(w) = |a(w)|2, a.{Lo) = j W{t){cy.{t))e'^'dt, (2a) 

{a{t)) = (V'(r,t)|Vy(r)+E(i)|V'(r,t)), (2b) 

where {a.{i)) is the dipole acceleration and S{uj) is pro- 
portional to the intensity of the emitted radiation at fre- 
quency u). Here W{t) is a standardly used temporal win- 
dow that prevents high-frequency artifacts at the bound- 
aries of the integration. 

For HJ, one observes a minimum in the spectrum of 
emitted radiation polarized in the x-direction because of 
interference between the two centers of the molecule Q . 
Using the plane-wave approximation for the returning 
electron, the (first and usually only observable) minimum 
occurs when the x-projection of the internuclear distance 
as seen by the returning wave packet is equal to half the 
de Broglie wavelength. Therefore the minimum will shift 
toward higher harmonics with higher angles between the 
laser polarization and molecular axis. Since it is a struc- 
tural minimum that depends only on the geometry of the 
bound state, one expects to see no shifts in the location 
of the minimum when different laser pulses are used or 
when instead Gaussian wave packets are used to generate 
harmonics in a laser-field-free electron- ion collision 

For the wave-packet simulation without laser field, the 
center of the potential is placed in the middle of the grid 
at (a;,y) = (0,0), and the Gaussian wave packet is in- 
troduced with its center at position (xq ~ 0), with 
the grid size denoted as y< Ly. The initial wave packet 
'0(r) is given by the superposition 

^(r) = V'o(r)+^'G(r), (3a) 

V>G(r) = Vcy^e"K^'(^-^)'+='''')+'^°^ (3b) 

where V'o(r) is the ground-state wave function, and c^, 
Cy quantify the momentum spread of the Gaussian wave 
packet in the x- and y-directions. The wave packet moves 



with a central momentum fcg < toward the molecular 
core. The norm of the Gaussian wave packet C should be 
set small to mimic the situation of HHG at the typically 
used intensities. We use C = 10~^. The momentum- 
spread parameter in the y-direction is chosen as 

Cy=rk\ko\, (4) 

where a tuning parameter is used to study the ef- 
fects of the different types of Gaussian wave packets 
and can be set to simulate the character of the contin- 
uum wave packet as generated by a laser pulse. The 
momentum-spread parameter in the x-direction Cx is set 
relatively large to allow for many harmonics to be probed 
by one wave packet. The propagation time is chosen such 
that a classical particle with momentum ko moves from 
{Lx/4:, 0) to {—Lx/4:, 0) during the propagation. As a re- 
sult, the strongest emission is expected at the middle of 
the propagation, such that little distortion is introduced 
when using a window function in the temporal Fourier 
transform for obtaining the power spectrum. As an ex- 
ample, for = 0.01 and ko = —1.78, the grid dimen- 
sions arc Lx ~ 383 a.u. and Ly — 1006 a.u. Here we use 
2304 X 6144 spatial grid points and 2000 time steps. The 
propagation time equals the optical period of a laser field 
with a 780-nm wavelength. 

For the simulation of the laser-induced HHG process 
we use a laser pulse linearly polarized in the x-direction. 
The initial state is set to the ground state tpoir). The 
time-dependent wave function is propagated for the laser 
pulse duration and two additional cycles after the end of 
the laser pulse to minimize distortions from the dipole 
acceleration window and to allow the wave packets to 
return to the nucleus. For the laser-induced calculations 
we use a grid measuring 280 x 84 a.u. with 1536 x 512 
grid points and 2000 time steps per optical cycle. 

III. HARMONIC INTENSITY 

The laser wavelength for the calculations of laser- 
induced HHG is 780 nm. For the Gaussian-wave-packet 
collisions, the same energy scale in units of harmonics 
of a 780 nm laser pulse is used. In Fig. [T] we show the 
alignment dependence of the emitted radiation for har- 
monic 49 polarized along the a;-direction. We compare a 
Gaussian wave packet with ko = —1.78 and = 0.01 to 
laser pulses with an intensity of / = 5 x 10^"* W/cm^ and 
different lengths. The laser pulses have a sin^ envelope 
of either three or five cycles, or a trapezoidal envelope of 
fifteen cycles length with five-cycle ramps. The carrier- 
envelope phase, i.e., the phase between the carrier wave 
and the envelope, is j for the shortest pulse and for 
the five-cycle sin^-pulse. The trapezoidal pulse has a si- 
nusoidal carrier wave. The intensity data are integrated 
over one harmonic order. The figure shows that the dif- 
ferent laser pulses give rise to minima that are near, but 
not exactly at, the same position as the minimum from 
the Gaussian wave packet. 
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FIG. 1. (Color online) Intensity of harmonic 49 vs alignment 
angle 6 for emission polarized along the x-direction. A Gaus- 
sian wave packet (black solid line) is compared to a three- 
cycle sin'^-pulse (red dashed line), a five-cycle sin^-pulse (blue 
dotted line) and a fifteen-cycle trapezoidal pulse (green dot- 
dashed line). 
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FIG. 2. (Color online) Location of the minimum in the 
alignment dependence of the intensity polarized along the x- 
direction. (Left) Blue crosses are for the Gaussian wave packet 
and red plusses for the three-cycle pulse of Fig. [1] (Right) 
Green squares are for the fifteen-cycle pulse of Fig. [1] and vi- 
olet triangles are for a ten-cycle trapezoidal pulse. The black 
solid line displays the two-cente r interference based on the 
SFA relation k{u}) = \/2{ui — Ip) and the black dashed line is 
the 7p-corrected result based on the relation fe(aj) = y/2Lu. 
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FIG. 3. (Color online) (Left) Same as Fig. [2] for a broad 
Gaussian wave packet starting far from the nucleus (blue 
crosses) and a broad Gaussian wave packet starting close to 
the nucleus (red plusses). (Right) A narrow Gaussian wave 
packet starting close to the nucleus (green squares) and a 
broad Gaussian wave packet starting far from the nucleus but 
with less energy (violet triangles). The blue crosses and red 
plusses have = 0.03 and ko — —1.78, the violet triangles 
have the same but kg — —1.38. The red plusses and green 
squares start at xq ~ 20 a.u., the blue crosses correspond to 
an initial position of xo = 95.8 a.u. and the green squares to 
xo ~ 74.1 a.u. 

of the minimum as a function of alignment angle 9. The 
laser pulses produce minima that follow the same trend 
as the Gaussian wave packet, but are scattered around 
the general trend. The results for short and longer pulses 
are scattered differently, but not less or more. The differ- 
ences between relatively long ten- and fifteen-cycle pulses 
pulses are small. This is expected, since both pulses 
are effectively almost cw-like. Even for the fifteen-cycle 
pulse the depletion of the ground state remains below 6%. 
The results suggest that when using the plane-wave ap- 
proximation for the returning electron in molecular imag- 
ing applications, a dispersion relationship in between the 
Lewenstein and /p-corrected relationships should be used 



A. Effect of the propagation step 



In Fig. [5] we plot the positions of the minima ^min in 
the alignment dependence versus harmonic order for dif- 
ferent laser pulses and the Gaussian wave packet from 
Fig. [TJ Also indicated in the figure are the curves that 
are predicted for the two-center minimum (R cos Omin = 
n/k) using either the energy-conserving relationship 
k{(jj) = y^2{u! — Ip) from the Lewenstein model [l^l, or 
using the /p-corrected relationship k{uj) ~ that has 
been adopted previously for molecular imaging [^. 
The physical argument for the /p-correction is that, when 
describing the returning electron as a plane wave, one 
should take into account that at the moment of recombi- 
nation its wave number is modified by the absorption of 
/p into the kinetic energy. From Fig. [5] it is clear that the 
Gaussian wave packet gives rise to a very smooth shift 



We compare the positions of interference minima for 
different types of Gaussian wave packets in Fig. [31 We 
vary the momentum spread of the wave packet in the per- 
pendicular direction, the central momentum of the wave 
packet and the position at which the wave packet starts. 
The striking and important observation is that all curves 
lie very close together. Apparently the position of the 
minimum is insensitive to the momentum distribution of 
the continuum wave packet. Only for the lowest harmon- 
ics we can observe some difference between the different 
kinds of continuum wave packets. Small differences ap- 
pear there between wave packets starting far away and 
those starting close to the nucleus, due to the effect of the 
long-range binding potential. Our findings suggest that 
in terms of the three-step model, the propagation step 
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has little effect on the observed position of the minimum 
and cannot account for the big fluctuations observed in 
Fig. EI 



B. Effect of the recombination step 

In the three-step model , the laser field during the re- 
combination step and possible interferences between dif- 
ferent parts of the continuum wave packet are ignored. 
In reality, the electronic wave packet does not recombine 
under laser-field-free conditions, and different trajecto- 
ries recombine at different times with different phases. 
To study the effect on the recombination process, we re- 
sort to a comparison of the minimum positions using ar- 
tificial pulses. The pulses are four-cycle sinusoidal pulses 
with a constant envelope corresponding to an intensity 
of / = 5 X 10"'^'* W/cm^, i.e., a section of a cw laser field. 
At t = 0, the electric field is £'(0) — 0. Optionally, we 
employ either or both of two methods to influence the 
recombination step: (i) setting the dipole acceleration to 
after some point in time during the propagation and 
(ii) turning off the laser fleld for the inner region near 
the nuclei after some point in time. 

Every half laser cycle, both a short and long classical 
electronic trajectory contribute to every harmonic peak 
(l3| . The distinction between short and long trajectories 
is based on whether the electron spends shorter or longer 
than 0.65r in the continuum, where T is the laser period. 
Setting the dipole acceleration to beyond t = Ta using 
a temporal width AT^ , Eq. (Pa|) becomes 



W'[t) = W{t)Sit), 

{1 for t < ti 

cos^(^j^^) f0lti<t<t2 

for t > t2, 



(5a) 
(5b) 

(5c) 



and 



where L is the propagation length, ti ~ — 

t2=Ta + ^. We use AT^ = O.IT. We set the dipole 
acceleration to at either Tq, = 0.95T, the return time 
of the most energetic trajectory, or at Tq, = 1.182r, the 
time at which the return momentum of the flrst half- 
cycle's long trajectory matches that of the second half- 
cycle's short trajectory. Thus with Ta = 0.95T we take 
into account only the short trajectories from the flrst half- 
cycle, and Tq, — 1.182T is the optimal point in time for 
selecting only a single pair of short and long trajectories. 

Additionally, we optionally turn off the laser in the 
inner region at time t = T\. To prevent artifacts, the field 
is turned off gradually in both space and time. The laser 
interaction is completely turned off for r — -^/x^ + < 
4, undisturbed for r > 6, and we use a sin^-transition 
between these two extremes. In the time domain, we use 
a smoothcned step function (convolution of a Gaussian 
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H = + V{x,y) + Z{x,y,t) ■ xE{t), 
Z{x, y, t) = Fix, y) + il- F{x, y))R{t), 

{0 for r < 4 

sin2(f(r-4)) for4<r<6 
1 for r > 6, 

R{t) = l(l- erf ({t - Ti) / (V2 • O.IT 



FIG. 4. (Color online) (Left) Same as Fig. [2]for a broad Gaus- 
sian wave packet starting far from the nucleus (blue crosses) 
and for harmonics generated from a sinusoidal laser pulse 
(red plusses). (Right) Harmonics from sinusoidal pulses with 
Tc = 1.182r (green circles) and with = 0.95T (violet 
points). 



with a step function) with a width of O.IT. In formula, 
the Hamiltonian in Eq. (jlbp is replaced by 



(6a) 
(6b) 

(6c) 
(6d) 



When Ti = 0.5r is used, this special setup allows us to 
compare near-physical harmonics to those generated in 
an identical setup where only trajectories starting dur- 
ing the first half-cycle contribute and with the laser fleld 
completely turned off during all recombinations. This 
method could be easily extended to filter out either the 
short or long trajectories. 

In Fig. 2] we compare the scattering around the 
Gaussian-wave-packet results from the sinusoidal pulse 
(red plusses) with those from setting the dipole accele- 
ration to at Tq, = 1.182T (green circles) and from set- 
ting the dipole acceleration to at Tq, = 0.95T (violet 
points). With Tq, = 1.182T the interference between the 
long and short trajectories leads to a strong, but reg- 
ular oscillation of the laser-induced results around the 
Gaussian- wave-packet results. The same interference be- 
tween the short and long trajectories can be seen in an 
associated harmonic spectrum as the top, black solid line 
in Fig.O When every harmonic peak is caused by a single 
trajectory (violet points in Fig. 2]), the interference dis- 
appears completely and the result is almost as smooth as 
that from the Gaussian wave packet. Under normal cir- 
cumstances, additional later returns from the same tra- 
jectories contribute to the spectrum. For a finite pulse 
length, different half-cycles also contribute differently be- 
cause of the pulse envelope. Additionally, for a numeri- 
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FIG. 5. (Color online) Harmonic intensity polarized in the 
x-direction for the green circles from Fig. |4] for 9 = 45° 
(black solid line) and 9 = 50° (red dashed line). Also for 
the violet points from Fig. [6] for 9 = 45° (blue dotted line) 
and 9 — 50° (green dot-dashed line). Solid and dashed 
arrows indicate spectral minimum positions predicted using 
k{uj) — \/2{uj — Ip) and k{uj) = V2uj, respectively. 
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FIG. 6. (Color online) Same as Fig. [2] for harmonics from a 
sinusoidal laser pulse with Ta = 1.182r (green circles) and 
additionally with Ti — 0.5T (violet points). 

cal calculation the dipole acceleration window W{t) also 
changes the contributions between different half-cycles. 
All of these together then smoothen but irregularize the 
oscillation of the green circles in Fig. 21 leading to the 
scattering of the laser-induced data points observed in 
Fig. [21 The interference between the short and long tra- 
jectories can probably also explain the scattering of the 
two-center minimum as a function of intensity as found 
by Gonoskov and Ryabikin [loj . 

The fact that the results for a single harmonic trajec- 
tory (violet points in Fig. [4]) lie so close to the Gaussian- 
wave-packet result, means that the laser field has no sig- 
nificant influence on the amplitude of the recombination 
matrix elements. This supports using HHG for molec- 
ular imaging 0, @, H, 0, as in a typical experimental 
setup only short trajectories contribute to the harmonic 
spectrum. Interestingly, however, in Fig. |6lwe show that 
turning off the laser field during the recombination docs 
have a significant effect on the interference between the 
short and long trajectories. In the figure, the green cir- 



cles are copied from Fig. |4l Additionally, we show the 
case where there is only a single set of short and long 
trajectories with additionally the laser pulse turned off 
in the inner region at Ti = 0.5T (violet points). The 
strong reduction in scattering amplitude for the violet 
points in Fig.jHlcan be understood from Fig. [51 where we 
plot the harmonic spectra for molecular alignment an- 
gles 45° and 50° for the case of the green circles (violet 
points) in Fig. |6l as the top (bottom) two curves. For 
the unmodified laser pulse we observe a significant shift 
of the trajectory interference positions in the harmonic 
spectrum when going from alignment at 45° to 50°. Al- 
though at first sight the bottom two curves in Fig.Olook 
more distinct from one another, a closer look reveals that 
the alignment dependence of the trajectory interference 
minima is actually a lot smaller with the laser field turned 
off in the inner region, as there is no shift visible. The 
strong scattering at low harmonics for the violet points 
in Fig. ini is caused by the fact that the finite temporal 
widths of the filters R{t) and S{t) suppress the complete 
lower end of the spectrum. 



IV. HARMONIC PHASE 

The two-center minimum in the harmonic spectrum is 
accompanied by a phase jump in the harmonic phase. Us- 
ing the plane- wave approximation, this should be a sharp 
TT-phase jump Q. However, in experiments a smaller and 
smoother phase jump is observed poj . Such deviations 
can be attributed to nonclassical momenta [2l| and to ef- 
fects of the Coulombic potential [12] . Similarly, a phase 
jump is observable when one considers a fixed harmonic 
as a function of 6. The phase of harmonic 49 for emis- 
sion polarized along x is shown in Fig. [71 The same set 
of laser pulses and Gaussian wave packet is used as in 
Fig. [1] The curves have been shifted such that for 6' = 
the phase is 0. The figure shows that both the Gaussian 
wave packet and the extremely short three-cycle laser 
pulse give rise to a mostly constant phase as a function 
of 6 with a phase jump slightly smaller than tt at the 
location of the minimum. The longer pulses have a more 
smeared-out phase jump. In the neighborhood of the 
minimum their jump is a lot smaller than tt but over 
the complete 6'-range the jump seems to be bigger than 
TT. This behavior for the longer pulses is probably an ef- 
fect of more, and longer, trajectories contributing to the 
harmonics. Different trajectories are associated with dif- 
ferent Coulomb corrections and therefore the harmonic 
phase becomes smeared out. This is in accordance with 
the shallower intensity minima in Fig. [1] for the longer 
pulses. 

It is interesting to investigate the phase jump for the 
different Gaussian wave packets of Fig. [31 This is plotted 
in Fig. ISl for a smaller range of 9 for clarity. The broad 
Gaussian wave packet starting closing to the nucleus (red 
dashed line) starts out with small perpendicular momen- 
tum components. Because of the short propagation time 
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FIG. 7. (Color online) Phase of the harmonic emission polar- 
ized along X for harmonic 49 generated by a Gaussian wave 
packet (black solid line) and generated by a three-cycle sin^- 
pulse (red dashed line), a five-cycle sin^-pulse (blue dotted 
line), and a fifteen-cycle trapezoidal pulse (green dot-dashed 
line). 
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FIG. 8. (Color online) Phase of the harmonic emission polar- 
ized along X for harmonic 49 generated by a broad Gaussian 
wave packet starting far from the nucleus (solid black line), a 
broad Gaussian wave packet starting close to the nucleus (red 
dashed line), a narrow Gaussian wave packet starting close to 
the nucleus (green dot-dashed line), and a broad Gaussian 
wave packet starting far from the nucleus but with less en- 
ergy (blue dotted line). 



before the interaction with the core, both Coulomb ef- 
fects on the momentum distribution and perpendicular 
momentum components will be relatively small in this 
scenario. As a result we observe a sharp almost- 7r-phase 
jump. To some extent, the difference with the narrow 
Gaussian wave packet (green dot-dashed line) is that 
the latter experiences roughly the same Coulomb effects, 
but starts out with much larger perpendicular momen- 
tum components. This leads to a much smoother phase 
jump. Independent of the central momentum of the wave 
packet, a broad wave packet starting far away from the 
nucleus (black solid and blue dotted lines) experiences a 
smoother phase jump because of Coulomb effects and the 
associated increased nonparallel momentum components. 



V. HARMONIC POLARIZATION 

High harmonic radiation is coherent with well-defined 
polarization [l^ . We can write the emitted radiation 
field Eom(i) as 



Eom (t) 



v/l + (e(a;))2 2tt 



(7) 

where c is a complex function and e is the ellipticity as 
a function of frequency ui. Here e can be either posi- 
tive or negative, e = represents linearly polarized light, 
and e = ±1 represent positively and negatively circularly 
polarized hght, depending on the choice of propagation 
direction. Furthermore, Bp is the unit vector in the main 
polarization direction, and is the unit vector in the 
perpendicular direction. The angle between the main 
polarization direction of the emitted radiation and the 
polarization axis of the laser pulse is (j). Similarly to the 
angle 6 between the molecular axis and the laser polar- 
ization axis, a counterclockwise rotation corresponds to 
a positive angle. However, we limit the definition of 
to the range [0,7r), whereas the complex number c cov- 
ers the full complex domain. At a given w, all possible 
types of coherent plane-wave radiation can be uniquely 
described by the parameters (0, e, c). 

Numerically the acceleration of the dipole moment is 
identified with the far-field harmonic field (up to an over- 
all phase). Therefore the Fourier transformed dipole 
acceleration is given by 



Oi{uj) 



(ep(w) + je(a;)eo(a;)) . (8) 



For an experimentalist the most practical way to mea- 
sure the polarization angle and ellipticity of the emitted 
radiation is to pass the harmonics through a polariza- 
tion filter and measuring the emitted intensity for many 
polarization directions. The direction of greatest emis- 
sion corresponds to the main polarization direction, and 
the emission in the orthogonal direction is a measure for 
the ellipticity of the emitted radiation [l^. Only under 
considerable effort for the experimentalist, the harmonic 
phase can be measured interferometrically as a function 
of alignment angle [l^, HI] or as a function of harmonic 



order [2J]. Theoretically one has easy access to the har- 
monic phase. Using the phase information, Eq. ([S]) can 
be inverted as described in the following to deduce the 
polarization parameters 0, e and the complex number c 
from the complex- valued and ay. As a measure of 
the ellipticity we define a cross term a as (we omit the 
w-dependence for simplicity) 



q;^; a,, sm 



(5 = 



(9) 



where 5 = arg(ay) — arg(aa;) and we used that |cp = 
jctp. The following equalities can be straightforwardly 
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FIG. 9. (Color online) Main polarization direction (j) for har- 
monic 49 generated by a Gaussian wave packet (black solid 
line) and generated by a three-cycle sin'^-pulse (red dashed 
line), a five-cycle sin'^-pulse (blue dotted line), and a fifteen- 
cycle trapezoidal pulse (green dot-dashed line). 



FIG. 10. (Color online) Main polarization direction (f) for 
harmonic 31 (black solid line), harmonic 41 (red dashed line), 
harmonic 51 (blue dotted line), and harmonic 61 (green dot- 
dashed line). The black crosses are the phase difference 5 for 
harmonic 31. A three-cycle sin^-pulse was used. 



derived for the polarization parameters in terms of a 
and (T, 



1-4 



tan ( 



a. 



y 



i eax 



ax + ieay' 

— -^Ux if la^;! > \ay\ 



else, 




(10a) 
(10b) 
(10c) 



sin e cos y 



where for c we picked the numerically most stable ex- 
pression. Numerically one will run into problems using 
the above conversion if the emitted radiation is either lin- 
early or circularly polarized. Therefore one should check 
beforehand if one of these conditions applies and use ap- 
propriate simplified conversion equations instead. 



A. Polarization direction 

The two-center interference minimum can also be ob- 
served in the polarization direction (p of the harmonics. 
Because the emission in the direction parallel to the laser 
polarization direction is strongly suppressed at the min- 
imum, we expect a -l-jump toward the minimum. The 
TT-jump for the harmonic phase in the a;-direction in Fig. [7] 
translates to a full 7r-rotation for cf). This is exactly what 
is observed for a Gaussian wave packet and harmonics 
generated by different laser pulses in Fig. [S] The polar- 
ization direction of the emitted radiation was averaged 
over one harmonic order using the total emitted intensi- 
ties as weights. 

The phase jump in the x-direction at the two-center in- 
terference minimum becomes smoother for low harmon- 
ics [20l - [2^ . In the following we investigate how the jump 
in the polarization direction depends on harmonic order. 



We plot the main polarization direction for different har- 
monics in Fig. [ini Here we plot the polarization data 
points corresponding to the exact harmonics, i.e., no av- 
eraging was done. The right-most curves for high har- 
monics show a simple polarization-direction jump around 
the two-center minimum. For the lowest harmonics (on 
the left-hand side) the behavior becomes more compli- 
cated. To explain this finding, we also plot the phase 
difference S between Uy and ax for harmonic 31 (black 

crosses). The intensity ratio |^ equals 1 for 9 = 16° 
and 9 = 37°, and it reaches a maximum of -pN- = 2 at 



9 ~ 25°. For a given ratio 



l"i<l 



the main polarization 



direction is aligned more along the laboratory-frame (x 
or y) direction with the higher amplitude if the phase dif- 
ference b between the laboratory-frame directions is far 
from or tt. We can observe this effect clearly in Fig.lTUl 
at 9 equal to 20°-25° for harmonic 31, where b is around 
^ and \ay \ is bigger than \ax\-, the relatively slow increase 
in the polarization direction shows the tendency that the 
polarization is clamped toward the y-direction (0 = ^). 

We fit the jump observed in Fig. [9] with a smoothened 
step function to determine the location 9^ and the width 
of the polarization-direction jump. When one plots 
the location 0p as a function of harmonic order, one 
obtains a result very similar to that shown in Fig. [2] 
In Fig. [TT] we plot the width A6' as a function of har- 
monic order for a Gaussian wave packet and different 
laser pulses. Again, the results for the different laser 
pulses are scattered around the Gaussian- wave-packet re- 
sult. The width of the jump for as a function of har- 
monic order does not depend on the laser pulse length. 
One observes that the width of the jump decreases with 
increasing harmonic order. An important part of this de- 
crease is due to a purely geometric effect: with increasing 
harmonic order, the minimum moves to higher 9 leading 
to a narrower interference pattern as a function of 9. 

In the absence of Coulomb effects, i.e., when ax and ay 
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30 40 50 60 70 80 

Harmonic order 

FIG. 11. (Color online) Width of the jump in (f). Blue crosses 
are for the Gaussian wave packet, red plusses for the three- 
cycle pulse, and green squares for the fifteen-cycle pulse. The 
brown line shows the curve expected from purely geometric 
considerations. 

have the same phases apart from 7r-jumps, the polariza- 
tion direction (j) is given by tan 4* ~ We assume that 
ay is ahgnment independent over the range of the two- 
center minimum, as also indicated by numerical tests that 
we have performed, and that the alignment dependence 
of ax comes purely from the two-center interference with 
the phase lag kRcos9, i.e., ax = ai°'' cos (ifci?cos 0) 
[2^ [25! . Then the variation A0 of the polarization direc- 
tion on varying 9 in the vicinity of 9p is proportional to 
A{kRcos9) ~ (Tr/cos^p) A(cos6'). Thus the width A9 
should be such that A(cos6') /cos6'p is independent of 

harmonic frequency, provided that Uy/a^x^ is frequency 
independent. The brown line in Fig.[TT]is obtained for A9 
if we set A(cos0) /cos^p arbitrarily equal to 0.185 using 
A (cos 6*) ~ cos {9p - ^) - cos {9p + ^) . A comparison 
between the brown line and the other curves shows that 
at the high end of the spectrum the decrease in A9 can- 
not be explained any more exclusively by the geometric 
effect. Since 9p varies very slowly in this range, we expect 

only a slow variation in ay/a^x"^ . This suggests that de- 
creasing Coulomb effects play a role, in accordance with 
Fig. [5J where for harmonics 50-80 we observe a transi- 
tion toward the curve predicted by the SPA dispersion 
relationship, also indicating decreasing Coulomb effects 
in this range. 



B. Ellipticity 

In Fig. [12] we plot the ellipticity e of the emitted ra- 
diation for harmonic 49 as a function of 9. Again the 
polarization data was averaged over one harmonic order 
using the intensities as weights. The plot shows that both 
a Gaussian wave packet and different laser pulses give 
rise to both significant and varying elliptical emission. 
Nonzero ellipticity means that the harmonics in the x 
and y-directions are emitted with different phases. Using 
the plane- wave approximation for the returning electron. 
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FIG. 12. (Color online) Ellipticity of the harmonic emis- 
sion for harmonic 49 generated by a Gaussian wave packet 
(black solid line) and generated by a three-cycle sin^-pulse 
(red dashed line), a five-cycle sin^-pulse (blue dotted line), 
and a fifteen-cycle trapezoidal pulse (green dot-dashed line). 
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FIG. 13. (Color online) Alignment angles of zero elliptic- 
ity. (Left) Blue crosses are for a broad Gaussian wave packet 
starting far from the nucleus and black triangles for a narrow 
Gaussian wave packet starting close to the nucleus. (Right) 
Red plusses for the three-cycle pulse and green squares for 
the fifteen-cycle pulse. 



one would not expect to see any ellipticity for a symmet- 
ric molecule (26j . Since the ionization and propagation 
step are identical for the two components of the radia- 
tion, the ellipticity must come from the recombination 
step. Therefore, this result confirms that the Coulomb 
effects can lead to significant ellipticity. The ellipticity 
for parallel or perpendicular alignment is zero, because 
at these alignment angles, the perpendicular component 
of the emitted radiation vanishes. The ellipticity goes 
through zero at some intermediate alignment angle. If 
we plot the angle of zero ellipticity as a function of har- 
monic order, we arrive at Fig. 1131 This plot shows that in 
the close vicinity of the two-center interference minimum, 
the ellipticity goes through zero. This is as expected, be- 
cause at the location of the minimum, the a;-component 
of the emitted radiation is very small. Because the x- 
component has opposite signs before and after the min- 
imum, the ellipticity changes handedness through the 
minimum. 
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FIG. 14. (Color online) Extrema of ellipticity as a function 
of harmonic order. Blue crosses represent a broad Gaussian 
wave packet starting far from the nucleus, black triangles a 
narrow Gaussian wave packet starting close to the nucleus, red 
plusses a three-cycle pulse, and green squares a fifteen-cycle 
pulse. 
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FIG. 15. (Color online) Alignment angles for extrema from 
Fig. [14l Same wave packets and pulses as in Fig. [14] 



For each harmonic, we can also plot the extrema of 
the ellipticity that can be reached and the alignment an- 
gles at which those extrema are reached. The results arc 
shown in Figs. [Til and [T51 respectively. We observe that 
for sufficiently high harmonic orders, the ellipticity ex- 
trema become smaller in absolute value and move closer 
to the two-center interference minimum with increasing 
harmonic order. Since the ellipticity is an indicator of 
non-plane- wave character, the decreasing ellipticity is an- 
other signature of decreasing Coulomb effects for higher 
harmonic orders, which correspond to higher return mo- 
menta. 



structural minimum from the two-center interference be- 
tween the two lobes of the orbital. The position of 
this minimum is reproduced using an effective plane- 
wave momentum that transitions from the Jp-corrected 



k{LL}) = v2w at low harmonics to the SFA-based fc(w) = 
y^2{uj — /p) at high harmonics. A laser-induced HHG 
spectrum shows the same behavior if only a single elec- 
tronic trajectory contributes per harmonic, as is the case 
for a typical experimental setup. This justifies using 
HHG for molecular imaging as the laser field has no sig- 
nificant effect on the amplitude of the recombination ma- 
trix element. When a single set of short and long electron 
trajectories contributes to the spectrum, the interference 
between the two trajectories causes a large but regular 
oscillation around the general trend. Introducing more 
and longer trajectories by using longer pulses has the 
effect of averaging out the oscillations to a smaller scat- 
tering around the Gaussian-wave-packet result and leads 
to a smoother interference minimum as a function of 9. 

Our results show that the effect of the Coulomb po- 
tential can lead to significant ellipticity of the emitted 
radiation. Around the interference minimum, the main 
polarization angle makes a 7r-jump and the ellipticity goes 
through zero. The Coulomb effects are less important at 
higher harmonics. Therefore we observe decreasing over- 
all ellipticity and a relatively sharp jump in the polariza- 
tion direction at the high end of the spectrum. 

In the wave-packet calculations, the Coulomb effects 
could be investigated in more detail by changing the 
strength of the potential for the evolution of the contin- 
uum part. This may be subject of future work. Finally, 
we mention that for randomly oriented molecules, the 
perpendicular harmonic components and thereby also the 
ellipticity vanish due to the cylindrical symmetry of the 
system around the laser polarization axis. The behavior 
of the phase is more complicated. Since, however, HHG is 
dominated in our case by the large orientation angles for 
geometrical reasons, we do not expect a clear signature 
of the phase jump for randomly oriented molecules. 
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